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Abstract. We study systems of globally coupled interval maps, where 
the identical individual maps have two expanding, fractional linear, onto 
branches, and where the coupling is introduced via a parameter - com- 
mon to all individual maps - that depends in an analytic way on the 
mean field of the system. We show: 1) For the range of coupling pa- 
rameters we consider, finite-size coupled systems always have a unique 
invariant probability density which is strictly positive and analytic, and 
all finite-size systems exhibit exponential decay of correlations. 2) For 
the same range of parameters, the self-consistent Perron-Frobenius op- 
erator which captures essential aspects of the corresponding infinite-size 
system (arising as the limit of the above when the system size tends to 
infinity), undergoes a supercritical pitchfork bifurcation from a unique 
stable equilibrium to the coexistence of two stable and one unstable 
equilibrium. 



1. Introduction 

Globally coupled maps are collections of individual discrete-time dynami- 
cal systems (their units) which act independently on their respective phase 
spaces, except for the influence (the coupling) of a common parameter that is 
updated, at each time step, as a function of the mean field of the whole sys- 
tem. Systems of this type have received some attention through the work 
of Kaneko [9l[T0] in the early 1990s, who studied systems of N quadratic 
maps acting on coordinates x\, . . . ,xn G [0,1], and coupled by a parameter 
depending in a simple way on x := iV~ 1 (a;i + • • • +xn). His key observation, 
for huge system size N, was the following: if (x )t=o,i,2,... denotes the time 
series of mean field values of the system started in a random configuration 
(x\, . . . ,a;jv), then, for many parameters of the quadratic map, and even for 
very small coupling strength, pairs (x t ,x t+1 ) of consecutive values of the 
field showed complicated functional dependencies plus some noise of order 
N^ 1 / 2 , whereas for uncoupled systems of the same size the x*, after a while, 
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are constant up to some noise of order N~ 1 / 2 . While the latter observation 
is not surprising for independent units, the complicated dependencies for 
weakly coupled systems, a phenomenon Kaneko termed violation of the law 
of large numbers, called for closer investigation. 

The rich bifurcation structure of the family of individual quadratic maps 
may offer some explanations, but since a mathematically rigorous investiga- 
tion of even a small number of coupled quadratic maps in the chaotic regime 
still is a formidable task, there seem to be no serious attempts to tackle this 
problem. 

A model which is mathematically much easier to treat is given by coupled 
tent maps. Indeed, for tent maps with slope larger than v2 and moderate 
coupling strength, a system of N mean field coupled units has an ergodic 
invariant probability density with exponentially decreasing correlations [13j . 
This is true for all N and for coupling strengths that can be chosen to be the 
same for all N . Nevertheless, Ershov and Potapov [7] showed numerically 
that (albeit on a much smaller length scale than in the case of coupled 
quadratic maps) also mean field coupled tent maps exhibit a violation of 
the law of large numbers in the aforementioned sense. They also provided 
a mathematical analysis which demonstrated that the discontinuities of the 
invariant density of a tent map are at the heart of the problem. Their 
analysis was not completely rigorous, however, as Chawanya and Morita [2] 
could show that there are indeed (exceptional) parameters of the system 
for which there is no violation of the law of large numbers - contrary to 
the predictions in [7]. On the other hand, references |17yi8j contain further 
simulation results on systems violating the law of large numbers. (But at 
present, a mathematically rigorous treatment of globally coupled tent maps 
that is capable of classifying and explaining the diverse dynamical effects 
that have been observed does not seem to be in sight either.) These studies 
were complemented by papers by Jiirvenpaa [8] and Keller [12] . showing 
(among other things) that globally coupled systems of smooth expanding 
circle maps do not display violation of the law of large numbers at small 
coupling strength, because their invariant densities are smooth. 

Given this state of knowledge, the present paper investigates specific sys- 
tems of globally coupled piecewise fractional linear maps on the interval 
X := [— 5,2]' wnere each individual map has a smooth invariant density. 
For small coupling strength, Theorem 4 in [T2] extends easily to this setup 
and proves the absence of a violation of the law of large numbers. For larger 
coupling strength, however, we are going to show that this phenomenon does 
occur in the following sense: 

Bifurcation: The nonlinear self- consistent Perron- Frobenius operator 
(PFO) P on Li(X, A), which describes the dynamics of the system 
in its thermodynamic limit, undergoes a supercritical pitchfork bi- 
furcation as the coupling strength increases. (Here and in the sequel 
A denotes Lebesgue measure.) 
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Mixing: At the same time, all corresponding finite-size systems have 
unique absolutely continuous invariant probability measures fiN on 
their iV-dimensional state space, and exhibit exponential decay of 
correlations under this measure. 

Stable behaviour: In the stable regime, i.e. for fixed small coupling 
strength below the bifurcation point of the infinite-size system, the 
measures hn converge weakly, as the system size N — > oo, to an 
infinite product measure (uq • A) N , where uq is the unique fixed point 
of P. 

Bistable behaviour: In the bistable regime, i.e. for fixed coupling 
strength above the bifurcation point of the infinite-size system, all 
possible weak limits of the measures are convex combinations of 
the three infinite product measures (u r • A) N , r G {— r*, 0, r*}, where 
now uq is the unique unstable fixed point of P and u± Tt are its two 
stable fixed points. (We conjecture that the measure (uq ■ A) N is not 
charged in the limit.) 

This scenario clearly bears some resemblance to the Curie- Weiss model from 
statistical mechanics and its dynamical variants. 

We also stress that a simple modification of our system leads to a variant 
where, instead of two stable fixed points, one stable two-cycle for P is cre- 
ated at the bifurcation point. This may be viewed as the simplest possible 
scenario for a violation of the law of large numbers in Kaneko's original sense. 

In the next section we describe our model in detail, and formulate the 
main results. Section [3] contains the proofs for finite-size systems. In Sec- 
tion [5] we start the investigation of the infinite-size system via the self- 
consistent PFO P. We observe that this operator preserves a class of prob- 
ability densities which can be characterised as derivatives of Herglotz-Pick- 
Nevanlinna functions. Integral representations of these functions reveal a 
hidden order structure, which is respected by the operator P, and allows us 
to describe the pitchfork bifurcation. In Section [5] this dynamical picture 
for P is extended to arbitrary densities. Finally, in Section [6j we discuss the 
situation when some noise is added to the dynamics. 



2. Model and main results 

2.1. The parametrised family of maps. Throughout, all measures are 
understood to be Borel, and we let P(B) := {probability measures on B}. 
Lebesgue measure will be denoted by A. We introduce a 1-parameter family 
of piecewise fractional- linear transformations T r on X := [—9,9], which will 
play the role of the local maps. To facilitate manipulation of such maps, we 
use their standard matrix representation, letting 

Jm(x) := — — — - for any real 2 x 2-matrix M = ( a \ ) , 
cx + a V c " / 
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so that f' M (x) = (ad — bc)/(cx + d) 2 and Jm ° In = Imn- Specifically, we 
consider the function Jm t , depending on a parameter r £ (—2,2), given by 
the coefficient matrix 



M r 



r + 4 r + 1 
2r 2 



One readily checks that fu r (-\) = ~\, Sm t {\) = f, /M r («r) = | for 
a r := — r/4, and that (the infimum being attained on <9X) 

/m 0) = o > 2 = inf f' M > forxGl. 

MA ' 2(rx + l) 2 2 + \r\ X Mr 

The latter shows that ju T is uniformly expanding if and only if |r| < |, and 
we define our single-site maps T r : X — > X with parameter r £ (—2/3, 2/3) 
by letting 

T r (x) := f Mr { X ) mod (z + ±) = { . _ . , . ™ ["K^' 

V 2/ [ lMr\ x ) - 1 — lNA x ) on (OV> 2 J ' 



where 



o i 1 



We thus obtain a family (7 1 r ) r . g (_ 2 /3 i 2/3) of uniformly expanding, piecewise 
invertible maps T r : X — ► X, each having two increasing covering branches. 
Note also that this family is symmetric in that 

(2.1) - T r (-x) = T_ r (x) for r £ (-|, |) and x £ X. 

According to well-known folklore results, each map T r , r £ (—2/3,2/3), 
has a unique invariant probability density u r £ 2? := {u £ A) : ti ^ 

0,J x udX = 1}, and T r is exact (hence ergodic) w.r.t. the corresponding 
invariant measure. Due to (|2.ip . we have u_ r (x) = u r (—x) mod A. We 
denote the Perron- Frobenius operator (PFO), w.r.t. Lebesgue measure A, 
of a map T by Pr, abbreviating P r := Pt t - In our construction below we 
will exploit the fact that 2-to-l fractional linear maps like T r in fact enable 
a fairly explicit analysis of their PFOs on a suitable class of densities. In 
particular, the u r are known explicitly: 

Remark 1. Let 7 r := j^p, d~ r := j^p. Then 

rS r i 2r 2 

(2.2) u r (x):= yz d V = -( T, ^7 n , ^ 

Jy r (I — xyY (rx — (1 — r))(rx — (1 + r)J 



is an integrable invariant density for T r , see |21j . Its normalised version 

/ r 2_ 4 \" 1 

(2.3) u r (x) := log • w r (x) 



is the unique T r -invariant probability density. 
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-0,-4 -0,2 0,0 0,2 0,4 00 -1 1 . i . 1 . i : r- 

x -0,4 -0,2 0,0 0,2 0,4 



Figure 1. The functions T r (left), and u T (right), for r = — \. 

The key point in the choice of this family of maps is that for r < 0, T r 
is steeper in the positive part of X than in its negative part, hence typical 
orbits spend more time on the negative part, which is confirmed by the 
invariant density (see Figure [1]). If r > 0, then T r favours the positive part. 

The heuristics of our construction is that for sufficiently strong coupling 
this effect of "polarisation" is reinforced and gives rise to bistable behaviour. 

2.2. The field and the coupling. For any probability measure Q G P(X), 
we denote its mean by 

(2.4) <f>(Q) := [ xdQ(x), 

Jx 

and call this the field of Q. With a slight abuse of notation we also write, 
for u G T>, 

(2.5) 4>{u) := f xu(x) dx = <p(Q) if u = 

Jx d\ 

and, for x G X , 

N N 

(2.6) 0(x):=-J>=0(Q) XQ = W Y, 6 *i- 

i=l i=l 

To define the system of globally coupled maps (both in the finite- and 
the infinite-size case) we will, at each step of the iteration, determine the 
actual parameter as a function of the present field. This is done by means 
of a feedback function G : X — > R := [—0.4, 0.4] which we always assume to 
be real-analytic Q and S-shaped in that it satisfies G'(x) > and G{— x) = 
—G(x) for all x G X, while G"{x) < if x > 0. The most important 
single parameter in our model is going to be B := G"(0) which quantifies 
the coupling strength. 



This is only required to obtain highest regularity of the invariant densities of the 
finite-size systems in Theorem [T] Everything else remains true if G is merely of class C 2 . 
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Remark 2. The following will be our standard example of a suitable feedback 
function G: 

(2.7) G(x) := Atanh (j[ x \ 

where < A ^ 0.4 and ^ B ^ 18. (This requires some numerical effort. 
For < A ^ 0.2 and ^ B ^ 15, elementary estimates suffice.) 

For the results to follow we shall impose a few additional constraints on 
the feedback function G, made precise in Assumptions I and II below. 

2.3. The finite-size systems. We consider a system T^ : X N — > X N of 
N coupled copies of the parametrised map, defined by (Tjv(x))j = T r r x \(xi) 
with r(x) := G((p(x.)). For the following theorem, which we prove in sec- 
tional we need the following assumption (satisfied by the example above): 

(2.8) Assumption I: G'{x) ^ 25 - 50|G(ai)| for all x G X. 

Theorem 1 (Ergodicity and mixing of finite-size systems). Suppose 
the S-shaped function G satisfies (|2.8p . Then, for any N G N, the map 
Tyv : X N — > X N has a unique absolutely continuous invariant probability 
measure [in- Its density is strictly positive and real analytic. The systems 
(T^r, /^tv) are exponentially mixing in various strong senses, in particular 
do Holder observables have exponentially decreasing correlations. 

The key to the proof is an estimate ensuring uniform expansion. After es- 
tablishing the latter in Section [3j the theorem follows from "folklore" results 
whose origins are not so easy to locate in the literature. In a C 2 -setting, 
existence, uniqueness and exactness of an invariant density were proved es- 
sentially by Krzyzewski and Szlenk [15]. Exponential mixing follows from 
the compactness of the transfer operator first observed by Ruelle [20] . For a 
result which applies in our situation and entails Theorem 1, we refer to the 
main theorem of |16j . 

2.4. The self-consistent PFO and the thermodynamic limit of the 
finite-size systems. Since the coupling we defined is of mean-field type, 
we can adapt from the probabilistic literature (see for example [22J2]) the 
classical method of taking the thermodynamic limit of our family of finite- 
size systems Tjv, as N — > oo. To do so, consider the set P(X) of Borel 
probability measures on X, equipped with the topology of weak convergence 
and the resulting Borel cr-algebra on P(X). Define T : P(X) — > P(X) by 

(2.9) T{Q):=QoT;^ Q) , where r(Q) := G(cf)(Q)). 

We can then represent the evolution of any finite-size system using T. 
Indeed, if £at(x) := ^Yli=i^xi is the empirical measure of x = (si) 1 ^ i ^ Ar , 
then 6 N : X N — > P(X) satisfies eyv T^r = T o e^. 
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Furthermore, when restricted to the set of probability measures absolutely 
continuous with respect to A, T is represented by the self- consistent Perron 
Frobenius operator, which is the nonlinear positive operator P defined as 

(2.10) P:Li(X,X)^Li(X,X), Pu := P G ^ {u)) u. 

Clearly, this map satisfies T(u ■ A) = (Pu) ■ A and preserves the set D of 
probability densities. Note, however, that it does not contract, i.e. there 
are u, v G V such that \\Pu — Pv\\l 1 (x,\) > \\ u ~ v llii(x,A)- 

One may finally join these two aspects, the action of T on means of Dirac 
masses, or on absolutely continuous measures, via the following observation: 

Proposition 1 (Propagation of chaos). Let Q = u • A € P(X), with u G 
T>. If (a^)i>i is chosen according to then, for any n 0, the empirical 

measures ejv(T^-(xi, . . . , xjv)) converge weakly to (P n u) ■ A as N — » oo. 

This result confirms the point of view that the self-consistent PFO P 
represents the infinite-size thermodynamic limit N — > oo of the finite-size 
systems Tjy- Its proof is reasonably simple (easier than for stochastic evolu- 
tions). The only difficulty is that T is not a continuous map on the whole of 
P(X). This can be overcome with the following lemma, which Proposition 
[I] is a direct consequence of, and whose proof is given in Section IA.1I 

Lemma 1 (Continuity of T at non-atomic measures). Assume that a 
sequence (Qn) n >\ in P(X) converges weakly to some non-atomic Q. Then 
(TQ n ) n>1 converges weakly to TQ. 

Here is an immediate consequence of this lemma that will be used below. 

Corollary 1. Assume that a sequence (vr n ) n>1 of T -invariant Borel proba- 
bility measures on P(X) converges weakly to some probability ir on P(X). If 
there is a Borel set A C P(X) with n(A) = 1 which only contains non-atomic 
measures, then tt is also T -invariant. 

2.5. The long-term behaviour of the infinite-size system. Our goal 
is to analyse the asymptotics of P on D. Some basic features of P can be 
understood considering the dynamics of 

#:(-§,§)-!*:= ^] , H(r) := G{4>{u r )), 

which governs the action of P on the densities u r introduced in § 12.11 as 

(2.11) Pu r = P H r r \U r . 

In studying P, we will always presuppose the following: 

(2.12) Assumption II: H is S-shaped. 

This assumption can be checked numerically for specific feedback functions 
G, like that of RemarkEl cf. gAjbelow. By pT]) . H(-r) = -H(r). Note, 
however, that r i— > 4>(u r ) itself is not S-shaped (see Figure [3]) so that the 
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S-shapedness of G alone is not sufficient for that of H. 

Assumption II will enter our arguments only via the following dichotomy 
which it entails: either 

H(r) has a unique fixed point at r = 
(the stable regime with H'(0) ^ 1 and r = stable), 

or 

H(r) has exactly three fixed points — r* < < r* 
(the bistable regime with H'(0) > 1 and ±r* stable). 

We will see that H' > and H'(0) = G"(0)/6, so that the stable regime 
corresponds to the condition G'(0) ^ 6. Observe now that 



(2.13) Pu r = u r iff 



r = (in the stable regime) 

r G {0, ir*} (in the bistable regime) 

(since u r ^ u r > for r ^ r', and each T r is ergodic). We are going to show 
that the fixed points uq = lx, and u± r * dominate the long-term behaviour 
of P on T> completely, and that they inherit the stability properties of the 
corresponding parameters — r* < < r*. Therefore, the stable/bistable 
terminology for H introduced above also provides an appropriate description 
of the asymptotic behaviour of P. 

Theorem 2 (Long-term behaviour of P on D ). Consider P : D — > T>, 
T> equipped with the metric inherited from Li(X,X). Assuming (I) and (II), 
we have the following: 

1) In the stable regime, uq is the unique fixed point of P, and attracts all 
densities, that is, 

lim P n u = uq for all u G T>. 

n— »oo 

2) In the bistable regime, {u~ rt , uq, n rt } are the only fixed points of P. Now 
Mo is unstable, while U- r - t and u rt are stable. More precisely: 

a ) u ±r, ore stable fixed points for P in the sense that their respective 
basins of attraction are L\-open. 

b) If u£ T> is not attracted by u^ r . t or u rt , then it is attracted by u$. 

c) uq is not stable. Indeed, uq can be L\- approximated by convex analytic 
densities from either basin. It is a hyperbolic fixed point of P in the 
sense made precise in Proposition^ of Section^ 

Example 1. In case G{x) = At&nh(Bx/A) with < A ^ 0.4 and ^ B ^ 
18, both theorems apply. The infinite-size system is stable iff B ^ 6, and 
bistable otherwise, while all finite-size systems have a unique a.c.i.m. in this 
parameter region. 

The theorem summarises the contents of Propositions [3j [4] and [5] of Sec- 
tion [5] (which, in fact, provide more detailed information). The proofs rest 
on the fact that PFOs of maps with full fractional-linear branches leave the 
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class of Herglotz-Pick-Nevanlinna functions invariant. This observation can 
be used to study the action of P in terms of an iterated function system on 
the interval [—2,2] with two fractional-linear branches and place dependent 
probabilities. In the bistable regime the system is of course not contractive, 
but it has strong monotonicity properties and special geometric features 
which allow to prove the theorem. 

Our third theorem, which is essentially a corollary to the previous ones, 
describes the passage from finite-size systems to the infinite-size system. 
Below, weak convergence of the fijy G P(X N ) to some /i£ P(X N ) means 
that f ipdfiN —* f fdfi for all continuous ip : — > K which only depend 
on finitely many coordinates. (So that f yd^N is defined, in the obvious 
fashion, for N large enough.) 

Theorem 3 (From finite to infinite size — the limit as N — > oo). 
The Tjv -invariant probability measures fxjsr of Theorem^ correspond to the 
T -invariant probability measures o e^ 1 on P(X). All weak accumulation 
points it of the latter sequence are T -invariant probability measures concen- 
trated on the set of measures absolutely continuous w.r.t. A. Furthermore: 

1) In the stable regime, the sequence (/ijy ° £ n 1 )n^i converges weakly to the 
point mass 5\. In other words, the sequence (hn)n^i converges weakly 
to the pure product measure A N on X N . 

2) In the bistable regime, each weak accumulation point it of the sequence 

(fiN ° % X )iV>i is °f the f° rm a &u- r *\ + (1 - 2q <W + a^ r ,A for some 
a 6 [0, |] . In other words, each weak accumulation point of the sequence 
(hn)n^i is of the form a(u_ rt \) N + (1 - 2a)\ N + a(u rt \) N . 

Remark 3. We cannot prove, so far, that a = ^, which is to be expected 

because uq is an unstable fixed point of P. In Section [6] we show that a = | 
indeed, if some small noise is added to the system. 

Proof of Theorem^ As ejv o Tjv = T o ejy, the T jv-invariant probability 
measures /xjy of Theorem [1] correspond to T-invariant probability measures 
UN ° e^ 1 on V(X). Their possible weak accumulation points are all concen- 
trated on sets of measures from V(X) with density w.r.t. A, see Theorem 3 
in [12] . (The proof of that part of the theorem we refer to does not rely on 
the continuity of the local maps that is assumed in that paper.) Therefore 
Corollary [1] shows that all these accumulation points are T-invariant prob- 
ability measures concentrated on measures with density w.r.t. A. In other 
words, they can be interpreted as P-invariant probability measures on T>. 
Now Theorem [2] implies that the sequence (/zjv ° e J^) N ^ converges weakly 
to the point mass 5 Uo \ in the stable regime, whereas, in the bistable regime, 
each such limit measure is of the form a S u _ r ^\ + (1 — 2a) 5 Uo \ + a 5 Urt \ for 
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some a G [0, ^] (observe the symmetry of the system). Now the correspond- 
ing assertions on the measures hn follow along known lines, for a reference 
see e.g. [12J Proposition 1]. □ 

3. Proofs: the finite-size systems 
We assume throughout this section that 

(3.1) \G(x)\ ^ 0.5 and G\x) ^ 25 - 50|G(x)| for all \x\ < -. 

In order to apply the main theorem of Mayer [16] we must check his as- 
sumptions (Al) - (A4) for the map T = T^. To that end define F : X N — > 

If by (P(x))i = /jif rW (xi). Obviously T(x) = F(x) mod (Z + %) N , 
and (Al) - (A4) follow readily from the following facts that we are going to 
prove: 

Lemma 2. F : X N — ► [— i, is a homeomorphism which extends to a 
diffeomorphism between open neighbourhoods of X N and [-j, %] N . 

Lemma 3. The inverse F _1 of F is real analytic and can be continued to 
a holomorphic mapping on a complex 5 -neighbourhood £1 of [— such 
that T _1 (f2) is contained in a 5' -neighbourhood of X N for some < 6' < 5. 

To verify these two lemmas we need the following uniform expansion 
estimate which we will prove at the end of this section. (Here |.| denotes 
the Euclidean norm.) 

Lemma 4 (Uniform expansion). There is a constant p G (0, 1) such that 
\\{D¥(yi))- 1 \\ < p for all N G N and x G X N . 

Proof of LemmalB Obviously F(X N ) C [— |, f ] ■ Hence it is sufficient to 
prove the assertions of the lemma for the map F := ^(F — (i, . . . , ^) T ) : 
X N — ► X N . As each fu r is differentiable on (—1,1) (recall that |r| < |), 
F extends to an analytic mapping from (—1,1)^ — ► K . By Lemma HI it is 
locally invertible on Q £ := (— ^ — e, ^ +e) for each sufficiently small e ^ 0. 
(Note that = int(X).) All we need to show is that this implies global 
invertibility of F\q q : Qq — ► £Iq, because then the possibility to extend F 
diffeomorphically to a small open neighbourhood of X N in WL N follows again 
from the local invertibility on f2 e for some e > 0. 

So we prove the global invertibility of F|^ : Qq — ► fio- As each fu r '■= 
|(/M r — h) '■ A — > X is a homeomorphism that leaves fixed the endpoints of 
the interval X, we have F(dX N ) C 8X N and F(Q ) Q flo- Observing the 
simple fact that Qq is a paracompact connected smooth manifold without 
boundary and with trivial fundamental group, we only need to show that 
F|n : is proper in order to deduce from [31 Corollary 1] that F|n 

is a diffeomorphism of Qq. So let A" be a compact subset of f2o- As F|q 
extends to the continuous map F : X N — > X N and as F(dX N ) C dX N , 
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the set F^iK) = F -1 (1T) C int(X N ) C X N is closed and hence compact. 
Therefore F|q : ^Iq — > is indeed proper. □ 

Proof of Lemma As F is real analytic on a real neighbourhood of X N , 
the real analyticity of F -1 on a real neighbourhood of [—5, follows from 
the real analytic inverse function theorem [14\ Theorem 18.1]. It extends to 
a holomorphic function on a complex 5- neighbour hood Q of [— ^ , |] - see 
e.g. the discussion of complexifications of real analytic maps in [144 pp. 162- 
163]. If 5 > is sufficiently small, Lemma H] implies that F _1 is a uniform 
contraction on $7. Hence the 5' in the statement of the lemma can be chosen 
strictly smaller than 8. □ 

Proof of Lemma^ Recall that (F(x))j = / r ( x )(xj) where r(x) = G(^(x)), 
0(x) = ^2 i= i Xi and we write f r instead of fu r - Denote g(x) = G'(^(x)), 
Ax(x) := diag^Ocx),...,/;^)) 

df r df r 

A 2 (x) := diag(— (xx), . . . , -^(x N )) 



/ N ■ ■ ■ N 



En 



\ N ■■■ N 

and observe that q(x) := (4 — r 2 )^-(x)/ f'.(x) simplifies to q(x) = 1 — Ax 2 
so that 

Ai(x) _1 A 2 (x) = - J 1 ^ diag (q(xi), q(x N )) =: - J L ^ 2 A 3 (x) . 
Then the derivative of the coupled map F(x) = (f r ( x )( x i)i . . . , f r (x)( x N)) is 
L>F(x) = Ai(x) + A 2 (x)E N g(x) = Ai(x) + A 3 (x)E N g(x) 

with 1 denoting the identity matrix. Letting 

q = q(x) := (q(xi), . . . ., q(x N )f 

A 1 vr 

e^v := (^••••'iV ) 

and observing that As(x) = diag(q(x)) so that /S.^{x)En = qe T w , the inverse 
of DF(x) is 

DF(x)- 1 = ( 1 - g(x) M X )e\) A^x)- 1 . 

V 4 - r 2 + 5 ((x)e T Ar q(x) J 

In order to check conditions under which F is uniformly expanding in all 
directions, it is sufficient to find conditions under which ||DF(x) _1 || < 1 
uniformly in x. Observe first that HA^x)" 1 !! sC (inf f' r )~ l < 1 for 

|r| < |. From now on we fix a point x and suppress it as an argument to all 
functions. Then, if v is any vector in W N , some scalar multiple of it can be 
decomposed in a unique way as av = q — p where p is perpendicular to q. 
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Denote p = e^p, q = e^q, and observe that q ^ as q has only nonnegative 

euclidian norm of (1 — - — H 

2 



entries. We estimate the euclidian norm of (1 — 4 _ T ,2_f„ p T n c l e Ar)( av ) : 



(3.2) 



1 9^ T Q^N I 

gp-gq \ 2 ,, ,, 2 2 



1 + , % r - q + p 

4 — r z + gq J 
-\ 2 

1 , II ||2 , || ||2 

1 + ^r— - q + p 

i + r- 1 ^ 2 ,,.,^ 2 



1 + 



q + p 



where T := 4 r . As p ^ iV 1 / 2 ||p|| and q ^ N 1 ||q|| 2 (observe that all 
entries of q are bounded by 1), we can continue the above estimate with 

^ II ,|2 m| „2^ 9 r-iAr-i/2,1 || r- 2 iv-i||q||2 

< q + p + 2 llp|lllq|l /, i t, i » r _i„ I,,,, + HpII 



(i + r- 1 iv- 1 ||q|| 2 ) 2 (l + r-^-iUqii 2 ) 2 

To estimate this expression we abbreviate temporarily t := iV~ 1//2 ||q||. Then 
^ t ^ 1, and straightforward maximisation yields: 

r-iJV-V2||q|| 9 



(l + r-iiv-iUqii 2 ) 2 ^ igv^Vt 
r-^-Miqii 2 1 



(l + r-iiv-iUqii 2 ) 2 ^ 4r 

So we can continue the above estimate by 



where ||q|| + ||p|| = I av I ■ Hence 



1 - r + ge T N q 
and therefore 



— qe T N ) (av) 

wq / 



1/2 

< ( 1+ 7=~1= "b TF 1 H av l 



1 2 + |r|\ / 9 1 \ 



1/2 



2 2- |r|y V i6>/3vT 4ry 

Observing T = 4 ~ r2 one finds numerically that the norm is bounded by 
0.99396 uniformly for all x, if —0.5 ^ r 0.5 and ^ g ^ 25 — 50r. Hence 
the map F is uniformly expanding in all directions provided (|3.1|) holds. □ 



w v( x ) ■= Ti 132 > xeX - 
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4. AN ITERATED FUNCTION SYSTEM REPRESENTATION FOR SMOOTH 

DENSITIES 

4.1. An invariant class of densities. The PFOs P r allow a detailed anal- 
ysis since their action on certain densities has a convenient explicit descrip- 
tion: Consider the family (w y ) ye (_2 2) of probability densities on X given 
by 

1 ~ 2/74 
(1 - xyf 

As pointed out in [2TJ (using different parametrisations), Perron- Frobenius 
operators Pt M of fractional-linear maps fu act on these densities via their 
duals fhf#, where 

m#: =(? J) M (° °) = (& «) ! °< M <:« 

in that 

( 4 -f) Pfu^j • w v) = (fj Wyd) ^) ' w y' with y ' = fM#{y)i 

for matrices M and intervals JCI for which /m(J) = X. (This can also 
be verified by direct calculation). Since //oi\(a;) = 4 , one can compute the 
duals 

0>U/) := f M *(y) = ~ f — TIT = 7 i -i \ — : —7 and 

/M r (£) {r + l)y + r + 4 
(4.2) y 

V ^ r ,v f f . 1 ^r(y) 2(y + r) 

My) ■= f N *(y) - 



f Nr (l) l-* r (v) (r-l)y-r + 4 

of the individual branches of T r , then express P r w y as the convex combina- 
tion 

P r w y = Pf Mr (l[_i, ar ) • + Pf Nr (l (ar ,i] • 
(4-3) =Pr(y) • ^(j,) + (1 -Pr(y)) • ^r r (y), 

with weights 

/Otr J _|_ y ^ ^ _|_ y 
dx = - - — and 1 - p r (y) = - + — 
-1/2 2 4 + ry 2 4 + ry 

It is straightforward to check that for every r £ (—2, 2) the functions <7 r , 
r r are continuous and strictly increasing on [—2, 2] with images <x r ([— 2, 2]) = 
[-2,2/3] andr r ([-2,2]) = [-2/3,2]. 

From this remark and (|4.3p it is clear that the P r preserve the class of 
those u £ V which are convex combinations u = f, 2 2 ^ w y dfj,(y) = j w,d/i 
of the special densities w y for some representing measure fj, from P(— 2, 2). 
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We find that P r acts on representing measures according to 

(4.5) P r I / w, dfj, = I w,d {C*n) with 

\J(-2,2) J J(-2,2) 

£*H := (p r ■ (J,) o a~ l + ((1 - p r ) ■ n) o r" 1 , 

where p r ■ fi denotes the measure with density p r w.r.t. [i. 

To continue, we need to collect several facts about the dual maps a r and 
r r . We have 

(4.6) a' r {y) = l 2(4 T r2) ^ 2 and r' r {y) - 



(ry + y + r + 4) 2 r (ry — y — r + 4) 2 ' 

showing that a r and r r are strictly concave, respectively convex, on [—2,2]. 
One next gets readily from (|4.2[) that l/cr r — l/r r = 1 wherever defined, and 

a r (y)<r r (y) for [_2,2]\{-r} 

while (and this observation will be crucial later on) a r and r r have a common 
zero z r := —r and 

, 2 2t 
(4.7) a r (z r ) = T T (z r ) = ^, and 0>(.z r +i) 



(r + l)£ + 4 



In the following, we restrict our parameter r to the set R = [— ^j, 
Direct calculation proves that, letting Y := [— |, |] , we have 

(4.8) a r {Y)UT r (Y) CY if r G -R, 

so that Y is an invariant set for all such a r and r r , and that 

(4.9) sup a' r = a' r (-^J ^ j and supr^ = T^Q^^^ for r £ R, 

which provides us with a common contraction rate on Y for the a r and r r 
from this parameter region. All these features of a r and r r are illustrated 
by Figure [21 

We denote by w(y) := 4>{w y ) the field of the density w y , and find by 
explicit integration that 

v ^ A IV l + y/2 1 A 1 (y\ 2k+1 

(4.10, = (j - ? ) log ^ + - = g (2fc+1)(2t + 3) (I) 

for y £ Y. In particular, uJ(0) = and wJ'(y) ^ g > 0, so the field depends 
monotonically on y. Note also that we have, for all fi £ P(Y), 

(4.11) (j) (J w,d/-^J = J J xwy(x) dx dfi(y) = J wdfi. 

We focus on densities u with representing measure [i supported on Y, i.e. 
on the class V := {u = f Y w,dfi: /i£ P(Y)}. Writing 

(4.12) r^-GUwd^j = G{<t>{u)) e R, 




Figure 2. The dual maps ay and r r for r = 0.3. The small 
invariant box has endpoints 7 r and r r . 



and recalling (|4.5p . we find that our nonlinear operator P acts on the rep- 
resenting measures from P(Y) via 



(4.13) P I J u/.d/xl = J w.d[C*nj with 

Z*n : = £ r * M n = (p rfl • /x) o a ~l + ((1 - p r J • /i) ° r" 1 . 

As supp(£*>), the support of is contained in oy(supp(/ii))U7y (supp(//)), 
it is immediate from (|4.8p that 

PV' C £>'. 

For r G i? = [— A, ^] we find that ay and r r each have a unique stable 
fixed point in Y, given by 

a r (^ r ) = 7 r := and T r [S r ) = 5 r : 



r + 1 r — 1 

respectively. Note that the interval Y r := [7r,^r] of width 2r 2 /(l — r 2 ) 
between these stable fixed points is invariant under both ay and the 
small boxed region in Figure [5J Furthermore, each of 7 r and 8 r is mapped 
to r under the branch not fixing it, i.e. 

oy(<5 r ) = r r (7 r ) = r, 

meaning that, restricted to Y r , ay and r r are the inverse branches of some 
2-to-l piecewise fractional linear map 

The explicit T r -invariant densities u r from (|2.3p can be represented as 
u r = f Y w, dfj, r with fj, r £ P{Y r ) C P(Y") given by 

(4.14) % (s)= f 1 „ g '- 2 -^" 1 



dX Kaj V 9r 2 -4y l-y 2 /4 
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Our goal in this section is to study the asymptotic behaviour of P onP', 
using its representation by means of the IFS C* . We will prove 

Proposition 2 ( Long-term behaviour of P on T>' ). Take any u € T>' , 
u = J Y Wy dfi(y) for some p € P(^)- The following is an exhaustive list of 
possibilities for the asymptotic behaviour of the sequence (C* n p) n ^o: 

(1) C* n p y 5q for some n ^ 0. Then (C* n p) n ^o converges to p, r „ and 
hence P n u converges to u rt in L\(X,X). 

(2) The interval conv(supp(£* n / u)) contains for all n ^ 0. Then 
(C* n p) n ^o converges to 5q and hence P n u converges to uo in Li(X, A). 
In this case also the length of conv(supp(£* n /x)) tends to 

(3) C* n p ~< 5o for some n ^ 0. Then {C* n p) n converges to p~ Tt and 
hence P n u converges to ti- r * in L\(X,X). 

In the stable regime, only scenario (2) is possible, so that we always have 
convergence of (C* n p) n -^o to Sq. 

Our arguments will rely on continuity and monotonicity properties of the 
IFS, that we detail below, before proving Proposition [2] in Section [4.41 



4.2. The IFS: continuity. Convergence in P(Y), limp n = p, will always 
mean weak convergence of measures, J Y p dp n — > f Y p dp for bounded con- 
tinuous p : Y — > R. Since Y is a bounded interval, this is equivalent to 
convergence in the Wasser stein-metric dm on P(Y). If and F v are the 
distribution functions of p and u, then 

/oo 
\F^(x)-F u (x)\dx. 
-oo 

The Kantorovich- Rubinstein theorem (e.g. [6j Ch.ll]) provides an additional 
characterisation: 

(4.16) dw(p,v) = sup / ifjd(p-v) 

</>:Lip y [i/>Kl JY 

for any p,u e P{Y). Here, Lip y [ip] := sup y:y , eY . y ^ y , \ip(y) - if){y')\ / \y - y'\ 
for any ifj : Y — ► R (and analogously for functions on other domains). We 
now see that there is a constant K > (the common Lipschitz bound for 
the functions y \— ► w y (x) on Y, where x £ X) such that 

(4.17) 



w,dp — I w,dv ^ K ■ dvp(/i, v). 

This means that, for densities from V , convergence of the representing 
measures implies Li-convergence of the densities. 
We will also use the following estimate. 

Lemma 5 (Continuity of (r,p) i— There are constants k%,K2 > 
such that 

dw(£lp, < Kidwiv, v) + K 2 \r - s\ 
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for all fi, v G P(Y) o,nd all r,s £ R. 

Proof. For Lipschitz functions tp on Y we have 

J ij)d{C* r ix - C*v) ^ Lip y [(ip o a r )p r + (tp o r r )(l - p r )) ■ d w (p,,u) 

(4.18) +sup Lip R [i>{o-.{y))p.{y) + tp(r.(y))(l - P.(y))] ■ V ~ s \ 

y eY 

Suppose that tp is C , with \tp'\ ^ 1. Then the first Lipschitz constant is 
bounded by 

°>) -PrlL + |K -Pr+Tr-(1 ~ Pr)^] < OO, 



Ki :-- 



sup [II (l 

r<Ei? 



and the second one by 



k 2 :-- 



sup 

r£R 



dr 



■Pr 



+ 



dr 



Pr 



+ 



dp r 

~ ar >'~dr~ 



< OO, 



and the lemma follows from the Kantorovich- Rubinstein theorem ()4.16p . 
since these C 1 functions ip uniformly approximate the Lipschitz functions 
appearing there. □ 

Corollary 2. The operators £*, r € R, and C* are uniformly Lipschitz- 
continuous on P(Y) for the Wasserstein metric. 

Proof. The Lipschitz-continuity of C* is immediate from Lemma [5j For C* , 
we recall from (|2.9p and f)4. 13j) that C*p, = C* fj, with r M = G(J Y wdfi) so 
that 

|r M - r„| ^ Lip(G) Lip(tTi) dw(/i,i/)> 



which allows to conclude with Lemma 



□ 



Remark 4. Rigorous numerical bounds give Ki ^ 0.5761 and &2 ^ 0.5334. 
These estimates can be used to show that not only the individual C* are 
uniformly contracting on P(Y), but also (under suitable restrictions on the 
function G) C* is a uniform contraction on P(Y*) where Y* is a suitable 
neighbourhood of the support of fi rr , and /j, rt is the representing measure 
of u rt with r* the unique positive fixed point of the equation r = G((p(u r )). 
Our treatment of C* , however, does not rely on these estimates, because it 
is based on monotonicity properties explained below. 

4.3. The IFS: monotonicity. On the space P(Y) of probability measures 
/i, v representing densities from D' , we introduce an order relation by defin- 
ing 

(4.19) fi ^ v :44> My £ Y : fi(y, oo) < u(y, oo) 

The symbols -<, ^ and y designate the usual variants of ^. We collect a 
few elementary facts on this order relation: 

(4.20) fj, ^ v if and only if f Y udfi J Y udv for each bounded and non- 
decreasing u : Y — > R. 
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In particular, if fj, ^ u, then J Y wdfj, ^ f Y wdi>, and hence ^ r v as well. 

(4.21) If n -< v and if pi,P2 : Y — > Y are non-decreasing and such that 
Pl(y) < P2(y) for all y, then p/J^i/o p~ 1 . 

(4.22) If X i/ and if J Y udp = § Y udv for some strictly increasing 
u : y — > R, then fj, = v. 

(4.23) Let z £ K. Then <5 2 ■< \x if and only if supp(/u) C [z,oo). 

We also observe that the representing measures fi r of the T r -invariant den- 
sities u r form a linearly ordered subset of P(Y). Routine calculations based 
on (|4.14p show that 

(4.24) fi r -< ii s if r < s. 

Our analysis of the asymptotic behaviour of the IFS will crucially depend 
on the fact that the operators C* and C* respect this order relation on P(V), 
as made precise in the next lemma: 

Lemma 6 (Monotonicity of (r, fj,) i— ► £*//). Let /i^g POO r,s £ R. 

a) If fi -< v, then C*fi -< £*is. 

b) If r < s, then C*/j, -< £*/i. 

c) If fj, -< v, then C*n -< C*v. 

Proof. Let u : Y — > R be bounded and non-decreasing and recall that 

(4.25) / ud(£*.fi) = / [u o o~ r ■ p r + u o r r ■ (1 — p r )\ dp,. 
Jy Jy 

[aj) Let fj, X v. In view of (|4.20p we can prove C*fi ^ L*v by showing that 
the integrand on the right-hand side of (|4.25p is non-decreasing. For this we 
use the facts that a r and r r are strictly increasing with a r < r r , and that p r 
is non-increasing, since p' r (y) = —jj^yi < 0. One gets then, for x<y, 

u(a r x)p r (x) + u{r r x){l — p r (x)) 
=u(o- r x)p r (y) + u(a r x) {p r (x) - p r (y)) +u{t t x){1 - p r {x)) 



(4.26) 



>o 

^U(a r x)p r (y) + u{T r x)(p r {x) -p r {y)) + u(t t x){1 - p r {xj) 

=u(a r x)p r {y) +u(t t x){\ -p r (y)) 
=u(a r y)p r (y) + u(r r y)(l - p r {y)) 

- [(u(o- r y) - u(o- T x))p r (y) + {u(r r y) - u{t t x)) (1 - p r (y))] 



>o 



Hence [i^v implies C* r [i -< C*v. Now, if u is strictly increasing, then (|4,26p 
always is a strict inequality, i.e. the integrand on the right-hand side of 
(|4.25j) is strictly increasing. Therefore, fj, -< v implies -< C*u by (|4.22p . 
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E]) We must show that (Q5l) is non-decreasing as a function of r. To this 
end note first that 

a 27) - y2 ~ 4 < o 

(4.28) 9 -^l= ( ,Q, 

(4.29) %M 8 - 2 " 2 A , 2 >0. 

ar (ry — y — r + 4) z 

Hence, if r < s, then 

u{a r x)p r {x) + u(r r x)(l — p r {x)) 
=u(a r x)p s (x) + ii(cr r x) (jp T {x) - p s {x)) +u(r r x)(l - p r (x)) 

" v ' 

>o 

^u(a s x)p s (x) + u(T s x)(p r (x) - p s (x)) + u(t s x)(1 -p r (x)) 
=u(a s x)p s (x) + u(t s x)(1 - p s {x)) , 
and for strictly increasing u we have indeed a strict inequality. 

[cj) This follows from[aj) and[b|): 

L* ii = Ct a -< £* v -< CZ v = L*v 

as r M = G(/ y wd/i) ^ G{J y wdv) = r v by (gJDJ. □ 

In §5.2l below. we will also make use of a more precise quantitative version 
of statement a). It is natural to state and prove it at this point. 

Lemma 7 (Quantifying the growth of [i ^ C*fi). Suppose that a,f3 > 
are such that v! ^ a and T' r ,a' r ^ (3. Then, for [i^v, 

(4.30) J ud{C* r v)- j ud(£* r n) > a/3 (J id du - J id d/A . 

Proof. Observing that u(a r y) — u(a r x) and u(r r y) — u(r r x) are ^ a(3(y — x), 
we find for the last expression in (|4.26p that 

[(u(a r y) - u(a r x))p r (y) + {u(r r y) - u{r r x)) (1 - p r (y))] > afiy ~ a(3x. 

This turns (14.260 into a chain of inequalities which shows that the function 
given by v(x) := u(a r x)p r (x) + u(r r x)(l — p r {x)) — aj3x is non-decreasing. 
Hence, by (14.201) . fj, ^ v entails J Y v d[i ^ j y vdv, which is (14.301) . □ 

4.4. Dynamics of the IFS and the asymptotics of P on V'. We are 

now going to clarify the asymptotic behaviour of C* on P(Y). In view of 
(|4.17j) . this also determines the asymptotics of P on D', and hence proves 
Proposition [21 

Our argument depends on monotonicity properties which we can exploit 
since the topology of weak convergence on P(V), conveniently given by the 
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Wasserstein metric, is consistent with the order relation introduced above. 
Indeed, one easily checks: 

(4.31) If (v n ) and (p n ) are weakly convergent sequences in P(Y) with 
Vn ^ v n for all n, then lim^ n ^ limP n . 

Recall from § 12.51 that, in the bistable regime, r* is the unique positive 
fixed point of the equation r = G(4>(u r )). For convenience, we now let 

:= in the stable regime. 

Then, in either case, u r with representing measure fj, r is fixed by P iff r £ 
{0,±r*}. By (|4.24|) we have /i_ r , ^ [1q = 5o ■< /i r „ with strict inequalities 
in the bistable regime. 

Lemma 8 (Convergence by monotonicity). If n ^ then [i -< 

£*fj, ^ £* 2 /x X . . . , and f/ie sequence (C* n fJ.) n ^o converges weakly to a mea- 
sure fj, r y fx with r £ {0, ±r*}. The same holds for y instead of ^. 

Proof. The monotonicity of the sequence (£*" \i) n ^a follows immediately 
from Lemma [Qbj) . Because of (|4.3ip . it implies that the sequence can have 
at most one weak accumulation point. Compactness of P(Y) and continuity 
of C* therefore ensure that (j£* n {j,) n ^o converges to a fixed point of £*, i.e. 
to one of the measures fi r with r £ {0, ±r*}, and (|4.3ip entails /x r >; /x. The 
proof for decreasing sequences is the same. □ 

The following lemma strengthens the previous one considerably. It pro- 
vides uniform control, in terms of the Wasserstein distance (|4,15p . on the 
asymptotics of large families of representing measures. 

Lemma 9 (Convergence by comparison). We have the following: 

a) In the stable regime, there exists a sequence (e n ) n ^o of positive real num- 
bers converging to zero such that 

d w {C* n fi, S ) ^ e n for n £ P(Y) and n £ N. 

b) In the bistable regime, for every y > there exists a sequence (e n )n^o °f 
positive real numbers converging to zero such that 

dyi/(>C* n /i, /i r _J < e n for [l £ P(Y) with /i y 5 y and n £ N. 

An analogous assertion holds for measures fj> ^ 8— y . 

Proof. As Y = [— |, |], we trivially have 5„ 2 /3 ^ A* ^ ^2/3 fo r an M £ P(^0- 
In particular, £*<5 2 /3 ^ <5 2/ / 3 , an d Lemma [5] ensures that (£* n £ 2 / 3 )n^o con - 
verges. Due to Lemma lObj) . we have 5q H /i r> = £* n / u r . is ^ C* n 5 2 /3 for all 
n ^ 0, showing, via (|4.3ip . that lim£* n 5 2 /3 = AV,- In the same way one 
proves that (£* n 5 -2/3)71^0 converges to //- r *- For the stable regime this 
means that both sequences converge to 5q = hq. 
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Hj) Assume we are in the stable regime. By the above discussion, 

E n := d H /(£* n 5_ 2 /3, S ) + <52/3i Sq) 

tends to zero. For any ll, (14.311) guarantees C* n S_ 2 /z ^ £* n /i H C* n S 2 / 3 for 
all n ^ 0. Hence Fc*n S _ 2/3 (y) ^ F^ nfi (y) ^ F^ n& ^{y) for all y, proving 

d W /(£* n ^,<5 ) < e n . 

[bj) Now consider the bistable regime. Note first that if there is a suitable 
sequence (e n )n>o for some y > 0, then it also works for all y' > y. Therefore, 
there is no loss of generality if we assume that y > is so small that 

(4.32) a G{w{y)) (y) = Q + ^) y + 0(y 2 ) > y 

(use (|4.10p . (|4.28p and (|4.2|> to see that this can be achieved.) Since 
C*S y = Pr{y)5 ar ( y ) + (1 - Pr(y))S Tr ( y ) and o>(y) < r r (y), we then have 
-< <^y -< ^G{w(y))^y = £*fiy> recan Lemma [6l Lemma [8] then implies that 
(vC* Tl 5 ?/ ) ri ^o converges to ^ r „. In view of the initial discussion, (C* n S 2 / 3 ) n ^o 
converges to /i r>fc as well, so that 

s n := cV(£* n <5y,AV*) + dw(>C* n 52/3)/"rJ 
defines a sequence of reals converging to zero. Now take any ll G POO with 
/i h S y , then £* n 5 2/3 ^ C* n Li h C* n 5 y for all n ^ 0, and dw(£* n //, /i r „) < e n 
follows as in the proof of[aj) above. □ 

This observation enables us to determine the asymptotics of C* n \i for any 
ll G POO which is completely supported on the positive half (0,2/3] of Y 
(meaning that ll y So, cf. (14.23f) ) . or on its negative half [—2/3,0). 

Corollary 3. Let ll G P(Y). 

a) In the stable regime, the sequence {C* n n) n ^o converges to So- 

b) In the bistable regime, if fi y So, then the sequence (£* n Lt) n ^o converges 
to Li rst . If ft -< So, it converges to /U_ r „. 

Proof, [aj) follows immediately from Lemma [Sfej) . We turn to[b]): Let r : = 
JyWoIli. Then r > because w > on (0,2/3] and ll y So- Therefore 
cr r (0) > <To(0) =0. Fix some y as in Lemma Iffi)]). w.l.o.g. y G (0, oy(0)). 
Then So -< 5 y ■< C*ll = C* ll since a r and r r map supp(/x) into [oy(0),2/3], 
so that indeed d\y(C* n Li, LL rt ) — ► as n — > oo by the lemma. □ 

It remains to investigate the convergence of sequences (£* n Li) n ^o when 
none of these measures can be compared (in the sense of -<) to So- To this 
end let [ao, bo] := Y. Given a sequence of parameters r\, r 2 , . . . G R define 

On := o- Tn o . . . o cr ri (ao) and b n := r r „ o . . . o r ri (6 ) 
for n ^ 1, and, for any ll = Lio G P(F) = P[ooi &o]j consider the measures 

r„ ° • • • ° L r x V- 
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Then supp(> n ) C supp(£^ n /i n _i) C a rn ([a n -i, 6„-i]) U r r „([a n _i, 6 n -i]) C 
[a n ,6 n ] by induction. Write [a,o] e := [a — e, 6 + e], where e ^ 0. The next 
lemma exploits the crucial observation that the two branches a r and r r have 
tangential contact at their common zero z r , see Figure [2j 

Lemma 10 (Support intervals close to zeroes). There exists some 
C G (0,oo) such that the following holds: Suppose that (r n ) n j>i is any given 
sequence in R. If for some £ ^ and n(e) ^ we have 

(# £ ) z rn+1 £ [a n ,b n ] e forn^n(e), 

then 

(4.33) lim \b n -a n \^Ce 2 , 

n— >oo 

3 

(4.34) lim max (\a n \ ,\b n \) ^ -e + Ce 2 , 

n— »oo 4 

and, in case e = 0, 

(4.35) G [a n) 6 n ] for n ^ n(0) + 1. 

Proof. Let e ^ and assume (Jfr e ). Note that, for n ^ n = n(e), 

if a n > z rn+1 , then < a„ + i < 3/4 • e, 
if b n < z rn+1 , then - 3/4 • e < b n+ i < 0, 
if a n ^ z r . n+1 ^ 6 n , then a n+i < < b n+ i. 

The first implication holds because = cr rn+1 (z rn+1 ) < cr r . n+1 (a 
oy n+1 increases strictly, and since by (Jl» £ ) we have a n G {z rn+1 ,z Tn+1 + s], 
whence a n+ \ < e-sup<T^ n+1 ^ 3e/4 due to (|4.9p . Analogously for the second 
implication. The third is immediate from monotonicity. 

Now, as a r and r r share a common zero z r , (|4.9p ensures b n+m — a n+m ^ 



\{bn+ m -i -an+m-i) hi case z r G [Ofl+m-i, &n+m-l]- Otherwise, note that z r 
is e-close to one of the endpoints, w.l.o.g. to a^+m-l- Since a r and r r are tan- 
gent at Zr, there is some C > s.t. < r r - +m (a n+m _i) - cr ra+m (on+ m -i) < 
^-e 2 in this case, while (|4.9p controls the rest of b n -^ m — a^_(- m . In view of 
diam(Y) = 4/3, we thus obtain, for m ^ 1, 

3 C 2 4 {3\ m 2 

"n+m ~~ a n+m ^ ^(On+m-1 ~~ a n+m-l) + ^ ^ 3 I 4 / + ^ £ ' 

Statement (|4.33p follows immediately. For the asymptotic estimate (|4.34p 
on max(|a n | , \b n \) = max(— a n ,b n ), use the above inequality plus the ob- 
servation that, by the first two implications stated in this proof, a n+m and 
—bn+ m never exceed 3e/4. Finally, if e = 0, (|4.35p is straightforward from 
(Jfr e ) and the third implication above. □ 

While the full strength of this lemma will only be required in the next 
subsection, the e = case enables us to now conclude the 
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Proof of Proposition^ The conclusions of (1) and (3) follow from Corol- 
lary [3j If neither of these two cases applies, then the assumption of (2) must 
be satisfied, and so condition (Jfro) °f Lemma [TUl is satisfied with n(0) = 0. 
Hence lim n ^ 00 max(|a n |, |fe n |) = by (|4.34j) . As the £* n u are supported in 
[a n ,b n ], these measures must converge to Jo- □ 

5. Proofs: the self-consistent PFO for the infinite-size system 

5.1. Shadowing densities and the asymptotics of P on D. We are 

now going to clarify the asymptotics of the self-consistent PFO on the set 
T> of all densities, proving 

Proposition 3 ( Long-term behaviour of P on V ). For every u G V, 
the sequence (P n u) n ^o converges in L±(X, A), and 

~ I = uo in the stable regime, 

lim P u < 

n-*co I £ {tf-r, j uo,u r:t } in the bistable regime. 

The basins {u G V : lim^^oo P n u = u± rt } of the stable fixed points u± rf are 
Li-open. 

(The set of densities attracted to no in the bistable regime will be dis- 
cussed in §5.21 below.) 

We begin with some notational preparations. Throughout, we fix some 
u £ T>. The iterates P n u define parameters r n := G((f)(P n ~ 1 u)) [n ^ 1). 
With this notation, P n u = P Tn . . . P ri u. 

We let 7Ttv, N ^ 1, denote the partition of X into monotonicity intervals 
of T rjv o . . . o T ri . Note that each branch of this map is a fractional linear 
bijection from a member of 7r/v onto X. Since the T r , r £ R, have a common 
uniform expansion rate, we see that diam(-7rjv) — > 0, and hence, by the 
standard martingale convergence theorem, M[u || o~(ttn)] — > u in Li(X, A), 
that is, 

(5.1) r] N := \\E[u \\ o-(ir N )] - u|| Li(x ,a) — >0 as N ~^ °°- 
Write 

4 N) : = P r N+k ■■■ p n (% II <r(n N )]) for k > and N ^ 1, 

and observe that ?;[ G V because it is a weighted sum of images of the 
constant function 1 under various fractional linear branches (recall (|4.ip and 

(3D). 

For N = we let Vq ^ := u r , and write, in analogy to the notation 
introduced for A r ^ 1, := P Th ...P ri (v^) and ?] := \\vq- 1 - u\\ Ll ( X ,\)- 
Obviously, vjj® G V for all k ^ 0. 
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Hence there are measures G POO sucn that = J Y w m dn^\ 

Observe also that 

(5.2) \\P N+k u - v[ N) \\ Ll{x ,x) < W for all k > and JV ^ 0, 
as ||P r || = 1 for all r, so that in particular 

(5.3) |#t,f>) - < r/iv, 1^(0(4^)) - r^ +fe+1 | < • ^ 
In addition, we need to understand the distances 

&k — llVf-fe - ^ U fc llii(A,A)- 
which, in fact, admit some control which is uniform in k: 

Lemma 11 (Shadowing control). There is a non- decreasing sequence 
(A n ) n ^o in (0,oo), not depending on u EV, such that 

(5.4) A^' k) < n N ■ A n for k,n^0 and N ^ 
Proof. Let r { n N,k) := G(</)(P n ~ 1 i;f ) )), and observe that (|3T7jl entails 



A lip p „( 7V ) p p „( 7V )| 



' 1 



Li(X,X) 



Applying Lemma [5] repeatedly, we therefore see that 
n-i 

K l |r/V+n+fc 



71-1 



r 

71 — 7 



7=0 
71-1 

= Ak 2 \G{cp{P N+n+k ^ 1 u)) - G(<f>(P n - i - 1 v^ f) ))\ 

7=0 

71-1 
7 = 

n-1 

^AUG'IU^^kI^ + A^), 

7 = 

where the last inequality uses (|5.2p . (Recall that P does not contract 
on Li(X,X), whence the need for the A^'^-term.) Letting K n := 1 + 

K IIG'Hoo K2 Y^i=o K \ an< ^ A^'^ := max.{Aj N ' k ^ : i = 0, . . . ,n — 1}, we thus 
obtain 

(5.5) A^ fc ) < K n ■ ( m + AS } ) 

which proves our assertion. □ 

We can now complete the 
Proof of Proposition 0. We begin with the easiest situation: 
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The stable regime. We have to show that lim n ^ 00 \\P n u — Uq\l,\(x,X\ = 0- 
Take any e > 0. Let (e n )n^o be the sequence provided by Lemma l9a| . and 
K the constant from (|4.17p . There is some n (henceforth fixed) for which 
Ke n < e/3. In view of (|5.ip . there is some iVo such that (1 + A„)r/yv < 2e/3 
whenever N ^ TVo- We then find, using (|5.2p . Lemma HTl and (|4.17p together 
with Lemma [QaD that 



\\P N+n u - u \\ Ll{XiX) < m + \\vW _ ?* V W \\ Li(XjX) + Ke n 

(5.6) < i lN + A n r] N + Ke n 

< e for N ^ N , 

which completes the proof in this case. 

The bistable regime. Given the sequence r n = G((j)(P n ~ 1 u)) as before, we 
let f 

a m b n ] C Y be the sequence of parameter intervals from Lemma 110] Ob- 
serve that the measures representing the Vn satisfy supp(//^) C [a n ,b n ] 
for all n and N. We now distinguish two cases: 

First case: For all e > we have from Lemma [TO) Then, for any 
e > 0, the lemma ensures that there is some n (henceforth fixed) with 
max (|a n | , \b n \) < e/AK, so that also dwil^n ,$o) < s/2K, whatever N. 
Due to (pn|) . r/Ar < e/2 for iV ^ iV , and we find, using (jOJ) and (|4"Tfj) , 

HP^^n - itolUi(X,A) < W + || - «o|Ili(X,a) 
< e for iV ^ iV , 

showing that indeed P n u — ► iio- 

Second case: there is some e > s.t. (Jfre) is violated in that, say, 

(5.7) z r „ < a n -i - e 

for infinitely many n. We show that this implies P n u — > u r . t . (If (4>e) is 
violated in the other direction, P n u — > ii_ rst then follows by symmetry.) 

In view of (|Q|) . and since (due to ^ <W+kJ (HSOD , and KW\ ) 

</>(. v k) ^ 4>{wa N+k ) = w(a N+k ), we have 

r w+fc+1 = G(<f>(P N+k u)) > G(0(vl N) ) - m ) 

^ G(<P(vl N) )) - WG'W^vn > G(w(a N+k )) - WG'W^n, 
and hence, observing that ||^r||oo ^ 1 and writing a{y) := o~G{w{y)){y) f° r 

(5.8) a N+k +i = o- rN+k+1 (aN+k) ^ &(a N+k ) - WG'W^m 
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for all N and k. Note that ef'(O) = a' (0) + ^^(0)^=0 • G"(0) • w'(0) = 
\ + \ ■ G'(0) • \ > 1, see (HID and (ITOD . Therefore, if we fix some 
U) E (l,cr'(0)), there exists some a* > such that cr(a) ua for all 
a £ (0, a*]. Without loss of generality, e/3 < a*. 

Now fix N such that HG'H^r/jv < (u> — l)e/3, and let iV + n + 1 satisfy 
E3). Due to dH7D, we have 



(5.9) iP'N+ri) > CTr N+n+1 (Zr N+n+1 + e) > e/3. 

Now, if ajy +n+ i ^ a*, then, by (15. 8j) . 

GJV+n+2 ^ <5"(a7V+n+l) — H^'Hoo 7 ^ 

^ or(o*) - (w - l)e/3 
^ wa* - (w - l)a* = a* > e/3. 
Otherwise, aAr +ri+ i € (0, a*), and again 

> we/3 - (u - 1) e/3 = e/3. 

It follows inductively that liminffca/j ^ e/3. More precisely: If N% and 
n\ are integers such that rj^ < g := [ | C?' [ 1 0CJ (u; — l)e/3 and ajvi+m > 
Zr Nl+ni+1 +e, then 

at > e/3 for k > N\ + n\. 

In particular, if the initial density u is such that t]q = \\u Tt — , u||l 1 (a',a) < f?> 
then we can take N\ = 0. 

Next, fix y := e/6 E (0,e/3), and choose a sequence (e n ) n ^o according to 
Lemma [£lb"|) . Then < y < and hence 5q -< 5 y ^ /li^ for A; > A?i + n\ so 
that the lemma implies 6w {£-* n l^^, Hr t ) ^ e n . Hence, by (|4.17p . 

Ijpn^W _ u ^ || Ll(X A ) ^ K ■ e n for > N\ + ni and all n, N. 
We then find, using (|5.2p and Lemma [TTT 

(5.10) ll^^-^lk^^^+ll^-W^lk^+^^n 

< r?Ar + A n 7] N + K -e n 

for > N\ + n\ and all n, AT. Now lim n _ >tX3 \\P n u — u rt \\li(x,\) = follows 
as in the stable case. 

It remains to prove that the basin of attraction of u r „ is Li-open. (Then, 
by symmetry, the same is true for u- rt .) As P is Li-continuous, it suffices to 
show that this basin contains an open Li-ball centered at u Tt , . To check the 
latter condition, first notice that z Tt < < supp(/i rt ) so that there is some 
ri\ > such that a^{ao) > 0. As we can assume w.l.o.g. that e < \z rt \, we 
have a™i(ao) > z Tx + e, and as P is Li-continuous, there is some g E (0, g) 
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such that a r = a rni o- • -oa ri (ao) > z r ni + £ whenever \\u — u r „\\ Ll ( X! \) < Q- 
Therefore we can continue to argue as in the previous paragraph (using the 
present n\ and N\ = 0) to conclude that linin^oo ||P n u — u rf ||l 1 (x,a) = 0- 

□ 

Remark 5. We just proved a bit more than what is claimed in Proposition [3j 
another look at equation (|5.10j) reveals that, in the bistable regime, the 
stable fixed point u rt of P is even Lyapunov-stable (and the same is true for 
U- Tt ). Indeed, fix e > 0, n\ G N and ~q > as in the preceding paragraph. 
That choice was completely independent of the particular initial densities 
investigated there, and the same is true of the choice of the constants K, A n 
and e n occuring in estimate (|5.10p . Now let 5 > 0. Choose n-i G N such that 
e„ 2 < and then 77 := min{g, 2 (i + a — ]}■ Then equation (|5.10p . applied 
with N = 0, shows that for each u G Li(X, A) with r/o = \\u— u ri , ||l 1 (x,a) < V 
and for each n ^ holds 

\\P ni+n2+n u - u r J\ Ll{x , x) ^ Vo (l + A n2 ) + K s n2 < 8 . 

5.2. The stable manifold of uo in the bistable regime. Let W s (uq) := 
{u G T> : P n u — > uq} denote the stable manifold of uq in the space of all 
probability densities on X. Clearly, all symmetric densities u (i.e. those 
satisfying u{— x) = u{x)) belong to W s (uq), because symmetric densities 
have field 4>{u) = so that also the parameter G{(p{u)) = 0, and symmetry 
is preserved under the operator Pq. 

However, W s {uq) is not a big set. In the present section we prove 

Proposition 4 (The basins of u± r ^ touch W s (uo) D T>'). Each density 
in W s {uq) n V belongs to the boundaries of the basins of u Tt and of u- Tt . 

We start by providing more information on the fields (fi(P n u) of orbits 
in W s (uq) fl V. Recall that for u = JyW.dn G V we have P n u = 
f Y w,d(£* n /i) (n ^ 0). Given such a density, we denote by R n (u) the 
"radius" of the support of £* n /j, i.e. R n (u) := inf{e > : supp(£* n /i) C 
[-£,£]}, and let <f> n (u) := <j>{P n u) = j Y wd(C* n n). 

Lemma 12 (Field versus support radius). In the bistable regime, for 
each u G W s [uq) n T>' there exists a constant C u > such that 

(5.11) \<f> n (u)\ ^ C u ■ (Rn(u)) 2 for n^0. 

Proof. In view of the explicit formula (|4.10p , we haveUJ'(O) = | and^"(0) = 
0, and therefore see that there is some e G (0, |) such that for every e G (0, e) 
and all y G [— 2e, 2e], 

(5-12) M < \w(y)\ < and \G(y)\ > (B - cs)\y\, 

where B := G'(0) > 6 and c, too, is a positive constant which only depends 
on the function G. In addition, elementary calculations based on (|4.2p and 
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(|4.6p show that letting k := max(l, ^p), £ can be chosen such that, for 
every e G (0, e) and r G [0, -Be), also 

K(2/)-^K& , \T r {y)--\ ^Be for |y| < e, 
(5.13) ^ ^ > <Tr{y) > Q - «el (y + r) ^ for y G [-r,e], and 
> r r (y) ^ a r (y) ^ (- + ej (y + r) > -Be for y G [-e,-r). 



(Recall that cr r and r r share a zero at z r = —r.) Finally, note that we can 

:*-*)(i+f-(* + §: 



w.l.o.g. take e so small that B := (§ - e) (l + § - (± + f)e) G (1,3]. (Due 



to Assumption I we have i? ^ 25. 

Consider some v = jyW» dv with v G P(^). We claim that for e G (0, -|) 
(5.14) |^(Pu)| ^ B ■ \4>(v)\ — e 2 if supp(i/) C [-£,e]. 

Denote r := G((p(v)) which by S-shapedness of G satisfies |r| < Be. In view 
of our system's symmetry, we may assume w.l.o.g. that r ^ 0. According 
to (|4TTll and (|4TT3|) we have 

<p{Pv) = (J) w, d(C*v)^j = j (w o ay) ■ p r dv + J (w o r r ) • (1 — p r ) dv 

so that, due to (joTH]) and (j5J3]) . 

(w o ay) • p r dv 

. /" ^ l[-B e ,o)°gr(y) Vbs] °<r r (y) \ , , / \ t / \ 

^ y y ( — ^r^2 — + 6 J • vrwMv) My) 

>i(^(i-) +i ^G-))- to+ ' ) ^ ) - 

Combining this with the parallel estimate for J y (uJ o r r ) • (1 — p r ) dv, we get 
Continuing, we find that 

I .Mrs -ydy{y)+ I ^—^■ y dv( y )+ 2 



Y 



2_ r) 6-6e 2 y^^ij 6 6 

2+£ ■ydv(y)+ j ^— ^ • y dv(y) + 
2m6 - be" 1 ./m 216 6 



j,o)°-oe- ^ 
|,o) ^[o, 3J 



> K-w{y)dv{y)+ K* ■ w(y) dv(y) + 

•/f-2,0) Jfo.il 



2 K£ 
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where K := (± + e)/(l-e 2 ) > K* := (± - «e)(l -e 2 ). As, because of (JEEJI , 
0(u) = J* y IZJ dz/ ^ | , so that r = G (4> (v)) ^ (B — ce) ■ 4> (v), we conclude 



0(P«) > f iT + ( ^ C£) ^) +(K-K*)J^ ^ w(y) du{y) 

1 \ / B fl c 



> <t>(v) i=:-e)[l + --[=- + -)e)-s 2 = B- 6(v) - e 2 , 



2 ) \ 6 \3 6, 
since K — K* $J 3e and ^ | whenever |y| ^ g ^ |. This proves (|5.14p . 

Now take any ti G W s (ito) nV . Then 4> n (u) — > 0, and the second alterna- 
tive of Proposition[2]applies, so that R n (u) ^ and (l+2Si2 n ,(n)) 2 ^ -^p- 
for all n larger than some 7%. In particular, 

(5.15) iWH 2 < (1 + 2Bi?„(u)) 2 i?„(u) 2 < ^±^R n (u) 2 

for these n in view of (|5.13p . Applying, for n n-g, the estimate (|5.14p to 
v := P n u and e := R n (u), we obtain 

|0 n+1 (tt)| ^ -B • \<j> n (u)\ - (R n (u)) 2 for n ^ rig. 

Suppose for a contradiction that (R n (u)) 2 < -^p- \4> n (u)\ for some n > nj. 
Then \<f> n+1 (u)\ > JB% n (u)|, and therefore (i2 n+ i(u)) 2 < ^{R n { u )) 2 < 

^rp-^Y^-\(fi n (u)\ < ^^-\(j) n+ i(u)\. We can thus continue inductively to see 
that |0 n (it)| < |0 n+ i(u)| < 1 4> n +2 (u) | < ... which contradicts (j) n (u) — > 0. 
Therefore \<t> n {u)\ ^ (R n {u)) 2 for all n > rig, and the assertion of our 
lemma follows. □ 

Lemma 13 (W s (uq) is a thin set for the order -<). In the bistable 
regime, if u = j Y w, du and v = j Y w, du are densities in V with u -< v, 
then at most one of u and v can belong to W s (uq). 

Proof. Suppose that u £ W s {uq). We are going to show that P n v — > u r „, 
i.e. C* n v — > jjL Tlt as n — > oo. 

Assume for a contradiction that also v 6 W s (uo). We denote the param- 
eters obtained from u by r n ^ := G((j)(P n ~ l u)) = G(J Y w d(£*( n ~^ u)) , and 
define r n>v analogously. Then our assumption implies that limn^oo r„ )M = 
lim^ — >OG 7*ri,v — 0* 

In view of (ETTD]) . w' ^ |, and one checks immediately that infy a' = || > 
| so that there is no > such that infy a' Tn ^ | for all n ^ no- Because 

of the strict monotonicity of C* (Lemma [6|) we have C* n ° \x -< C* n °v, so that 
(replacing /i and f by these iterates) we can assume w.l.o.g. that no = 0. 
Denote C*^ 71 ' := C* n o ■ ■ ■ o £^ ^ so that jC*™^ = C^i^fJ, and (/i r M being 
non-decreasing) £* n ^ ^ v for n ^ 1. Therefore 



30 JEAN-BAPTISTE BARDET, GERHARD KELLER AND ROLAND ZWEIMULLER 



> G Uwd{C^v)\ - G ^wd(C< n ^) 

> inf G> ■ ^wd{C<^v) - j^wd(C;^^ 



In view of the lower bounds for w' and a' Tn , , repeated application of 
the estimate (|4.30p from Lemma [7] yields 

inf G 1 I 

(5-16) r n ^ - r nijU ^ - X gn id d{v - p). 

Observe that the last integral is strictly positive because /j, -< v, cf. (14. 22ft . 
On the other hand, due to Proposition [2] there are e n \ such that 

supp(£*» U supp(£*V) C [-e n ,e n ], 

and as a' (0) = \ < § and r n ^r n>lJ -> (whence also z rn ^,z rni/ -> z = 0), 
there exists a constant C > such that e n ^ C'(|) n for n ^ n'. Hence 
| </> n (it) |, |0„(u)| < max{C u ,C„} • C 2 (|y) n for n ^ n! by Lemma [T2| and as 
r n ,u — r n ,fi ^ sup : uJ / -(|0 n (u)| + |(^ n (?;)|), this contradicts the previous estimate 

iprm]) . ' □ 

We can now conclude this section with the 

Proof of Proposition \4\ Suppose that u = j w,dfi £ W s {uq). For t £ (0, 1) 
let := j y w.cZ((l - t)fj, + tS 2 / 3 ) G £>'. Then !>- u, hence u® 
VF s («o) by the previous proposition. Therefore, due to Proposition [2] and 
monotonicity of £*, for any t, P n u^ converges to u Tt y uq as n — > oo. 

On the other hand, lim^o ||u — ||l 1 (a',a) = 0, so « is in the boundary 
of the basin of u rf . Replacing <5 2 /3 by <5_ 2/ / 3 yields the corresponding result 
for the basin of tt- r , ■ D 

5.3. Differentiability of P at C 2 -densities. As P is based on a parametrised 
family of PFOs where the branches of the underlying map (and not only their 
weights) depend on the parameter, it is nowhere differ entiable, neither as 
an operator on L±(X, A) nor as an operator on the space BV(X) of (much 
more regular) functions of bounded variation on X. On the other hand, as 
the branches of the map and their parametric dependence are analytic, one 
can show that P is differentiable as an operator on the space of functions 
that can be extended holomorphically to some complex neighbourhood of 
ice. 

Here we will focus on a more general but slightly weaker differentiability 
statement. 

Lemma 14 (Differentiability of P at C 2 -densities). Let u & C 2 (X) be 
a probability density w.r.t. A and let g £ Li(X,X) have J x gdX = 0. Then 

(5.17) Ap(u + rg)\ T=0 = P r (g) + w r (u) ■ G'{<t>{u)) <j>{g) 
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where r = G(cj>(u)), w r (u) := P r ((uv r )'), and v r (x) = 4 £_~>2~ ■ If we consider 

P as an operator from BV(X) to Li(X,X), then P is even differ entiable at 
each probability density u G C 2 (X) C BV(X) and 

(5.18) DP\ U = P r + G'((f)(u))w r (u)®(j) . 

Proof. In order to simplify the notation define a kind of transfer operator 
L by Lu := u + u o /a 1 ) and note that (Lu)' = Lu'. Observing that 

/at- 1 = /a/" 1 ° /r J1 V we have P r u = L ( u ° /m" 1 " /lf-0- Defme 

r r \0 1/ r lvl r 

( B \ At 2 — 1 

(5-19) v r (x):= -f M A(f Mr ( x )) 



Qj.J^J VMrWJ ~ 4 _ r 2 • 

For a function u G C 2 (X) denote by U the antiderivative of u. Then 
(5.20) 



where 



R s , r {x) := J\ s - t )^(U(f M -t(x)))dt . 
As |r(C/ o / M _i) = u o • |:/ M -i = {uv r ) o we have 

Together with (j5T2"0|) this yields 



P s u -P r u = L[uo f Mrl • /^.i -no f M -! • /^_! 

= (s - r) L ([uv r y o f M ^ ■ f M _^ + L 
= {s-r)P r ((uv r y)+L(R' Sir ) 

and |L(i?^ r )(x)| ^ C (s — r) 2 with a constant that involves only the first 
two derivatives of u. 

Now let u G C 2 (X) be a probability density, and let # G L\{X, A) be such 
that J gdX = 0. Let r := G{(j>{u)) and s := G(<£(tt + 5)). Then 

2 P(u + 9) - = (^.u - Pru) + P r g + (P s g - P r g) 

= (s-r) P r ((uv r )') + P r (g) + (P s g - P r g) + L{R' s r ) . 

This implies at once formula (|5.17p for the directional derivative, and as 
\\P s g — P r g\\i — ► (s — > r) uniformly for g in the unit ball of BV(X), also 
(|5.18p follows at once. □ 
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Proposition 5 ( u = 1 is a hyperbolic fixed point of P). In the bistable 
regime, u = \ is a hyperbolic fixed point of P\dc(b\{x) ^ n the following sense: 
the derivative of P : PnBV(A) — ► L\(X, A) at u = 1 has a one- dimensional 
unstable subspace and a codimension 1 stable subspace. 

Proof Let Q := DP\ u =x. As G"(0) = B and io (l) = P [2x] = [x], it follows 
from (15. 18ft that Q = Pq + B [x]®(f). (Here [2x] denotes the function x t— ► 2x, 
etc.) Observe now that = ^. ThenQfx] = P N + ^N = (| + ^)N 

so that, for B > 6, Q has the unstable eigendirection [x] with eigenvalue 
A := \ + y| > 1. On the other hand, as (f)(1) = 0, we have Ql = PqI = 1, 
so the constant density 1 is a neutral eigendirection, and finally, for / G 
ker(0) n ker(A), we have Qf = P f, so Var(Q/) < \ Var(/). □ 

6. The noisy system 

In Theorem [3] we proved that, in the bistable regime, each weak accumu- 
lation point of the sequence (/xjv £Ji~)n^i is of the form a6 u _ rif \ + (1 — 
2a)5 Uo \ + a5 Urt \ for some a G [0, |], i.e. that the stationary states of the 
finite-size systems approach a mixture of the stationary states of the infinite- 
size system. It is natural to expect that actually a = i, meaning that any 
limit state thus obtained is a mixture of stable stationary states of P. While 
we could not prove this for the model discussed so far, we now argue that 
this conjecture can be verified if we add some noise to the systems. 

At each step of the dynamics we perturb the parameter of the single-site 
maps by a small amount. To make this idea more precise, let 

r(Q, t) = G((j)(Q) + t) for Q G P(X) and t £ R , in particular 
(6-1) M 

r(x, t) = G(0(x) + t) for x G X and tel. 

Let rji, rj2, ■ ■ ■ be i.i.d. symmetric real valued random variables with common 
distribution g and \rj n \ ^ e. For n = 1,2,... and x G X N let us define the 
A^-valued Markov process (£n) n eN by £o = x and 

( 6 - 2 ) £n+l = T r(£„,»7«+i)(&0- 

Assume now that the distribution of £ n has density h n w.r.t. Lebesgue 
measure on X N . Then routine calculations show that the distribution of 
£ n+ i has density f R PN,th n dg{t) where P^,t is the PFO of the map Tjv,t : 
X N — > X , (Tjv,t(x))j = r r ( xt )(sci). It is straightforward to check that, 
for sufficiently small e, Lemmas [2] - U] from Section [3] carry over to all TVt 
(|t| ^ e) with uniform bounds, and that j^f, \PN,tf —PN,of\d\ N ^ constAr-e- 
Var(/) so that the perturbation theorem of [11] guarantees that the process 
(£n)neN has a unique stationary probability ^n,e whose density w.r.t. \ N 
tends, in Li(X N , \ N ), to the unique invariant density of Tn as e — > 0. This 
convergence is not uniform in N, however. Nevertheless, folklore arguments 
show that there is some e > such that, for all e G (0, e) and all iV G N 
the absolutely continuous stationary measure hn,£ is unique so that the 
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symmetry properties of the maps T r and the random variables rj n guarantee 
that Hn,e is symmetric in the sense that its density h,N,e satisfies h^^ e (x) = 

hN,e(-x). 

On the other hand, for each fixed e > 0, all weak limit points of the 
measures Hn,e ° e^y 1 as N — > oo are stationary probabilities for the P{X)- 
valued Markov process (H n ) nG N defined by 

compare the definition of T : P(X) — > P(X) in (|2,9p . The proof is completely 
analogous to the corresponding one for the unperturbed case (see Lemma Q] 
and Corollary [1]) . For e G (0, i) the symmetry of the Hn,e carries over to 
these limit measures Q in the sense that Q(A) = Q{fi : fi € A} for each 
Borel measurable set A C P(X) where £i(U) := fj,(—U) for all Borel subsets 
U C X. 

The following proposition then shows that, in the bistable regime and for 
small e > and large N, the measures hn,e are weakly close to the mixture 
| ((u_ r>t A) N + (u r:t A) N ) of the stable states for P; compare also Theorem 

Proposition 6 (Invariant measures for infinite-size noisy systems). 

Suppose G'(0) > 6 so that we are in the bistable regime and recall that the 
r] n are symmetric random variables. 

Then, for every 5 > there is eq > such that for each e E (0,£o) the 
stationary distribution Q e of H n on P(X) is supported on the set of mea- 
sures u ■ A G P(X) which have density u = j Y w, dfi £ D' with representing 
measures fi G P(^) satisfying dvy(/i, ^(/^-r* + /U r J) ^ S. 

Sketch of the proof. Let Q be a stationary distribution of H n that occurs as 
a weak limit of the measures Hn,e- So Q is symmetric. Just as in the proof 
of Theorem [3j where the "zero noise limit" , namely the transformation T 
is treated, one argues that Q is supported by the set of measures u ■ A, 
u G T). Arguing as in the derivation of (|5.2p one shows that densities in 
the support of Q can be approximated in L\(X, A) by densities from T>', 
and the stationarity of Q implies that Q is indeed supported by measures 
with densities from T>'. Therefore the process (S n ) n ^o can be described by 
the transfer operator £* of an iterated function system on Y just as the 
self-consistent PFO P is described by the operator C* in equation (I4.13|) . 
The only difference is that in this case one first chooses the parameter r ran- 
domly, r = G(J Y w dfj, + r/ n+ i) and then the branch o r or r r with respective 
probabilities p r (y) and (1 — p r (y))- 

Let y > be such that F{r/ n > y} > 0. Suppose now that for some 
realisation of the process (H n ) n ^o the numbers r(E n ,r) n +i) satisfy condition 
(J^e) °f Lemma [TU] for all e > 0. Then it follows, as in the proof of Proposi- 
tion [2J that lmin^oo r n (E n , rj n+ i) = and the measures H n converge weakly 
to A so that also limn^oo r(H n , 0) = 0. As r/ n > y > for infinitely many n 
almost surely, both limit cannot be zero at the same time, and we conclude 
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that almost surely there is some e > such that (Jfr £ ) is not satisfied. In 
particular, there are e > and n G N such that (4>e) is violated for n = n — 1 
with some positive probability k. 

Let E n = h n ■ A with h n = j Y w, dv n . (So h n and v n are random objects.) 
As in (15. 9ft we conclude that supsupp^) < — e/3 or inf supp(z/jj) > e/3 in 
this case. Without loss of generality we assume that the latter happens with 
probability at least t|. 

Next, as in (|4.32|) we may choose y £ (0, e/3) so small that < y < 
Vi := a G(w(y))(y) ^ inf supp(£*<5 2/ ). Hence, for reasons of continuity, there is 
£i > such that also y ^ inf supp(C*5 y ) if e G [0,£i). Therefore, in view of 
the monotonicity of the operator £*, we can conclude that inf supp(f n ) ^ y 
for all n ^ n with probability at least ^. Now fix 5 > 0. By Lemma 00 there 
is some (non-random) n± £ N such that dv(/(>C* ni f n , /i r ,) ^ | for all n ^ n 
with probability at least t|. But then, by continuity reasons again, there is 
£o £ (0, £i) such that dvi/(z/«+ni! /z r „) < S for all n ^ n with probability at 
least ?|. The claim of the proposition follows now, because (H n ) n is a Markov 
process and because the stationary distribution Q is symmetric. □ 



Appendix A. Some technical and numerical results 

A.l. Proof of Lemma HI It suffices to prove the convergence for evalua- 
tions of any Lipschitz continuous function <p defined on X. Let us denote 
fn — f{Qn) (resp. r = r(Q)), and a n (resp. a) the discontinuity point of 
T Tn (resp. T r ). Recall that a n = — (resp. a = — |). 

Let us fix £ > 0. Q being non-atomic, there exists 5 > such that the in- 
terval U := [c— (5, c+J] is of Q-measure smaller that £. The weak convergence 
of Q n to Q implies that c n tends to c, and that limsup n ^ +00 Q n (U) ^ Q(U). 



Let us choose no such that for all n ^ n , 
then has 



<pd{TQ)- / <MTQ n ) 

A 



|c n — c| < | and Q n (U) < e. One 



ipoT r dQ - / 99 o T rn dQ„ 

A Ja 



(A.l) 



/ (foT r d(Q - Q n ) + / (<p o T rn - Lp o T r ) dQ n 
Jx Ju c 

+ / (<poT rn - LpoT r )dQ n 
Ju 

I ipoT r d(Q-Q n ) + Lip(v?) sup |T rri — T r | + 2£||( / 9|| 0< 



Since the application cp o T r as a single discontinuity point, which is of zero 
Q-measure, the first term converges to zero. The second one also goes to 
zero since it measures the dependence of T r on its parameter away from the 
discontinuity point (one can make an explicit computation). 



A. 2. The fields of the densities u r . We start with some observations on 
the function ip{r) := 4>{u r ) that are based on symbolic computations and on 
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Figure 3. The functions ip(r) := <j){u r ) (left), Q(r) := 
(centre), and H(r) = Atanh(^d>(u r )) with A = 0.4 
and B = 8 (right). 



numerical evaluations. One finds 



lo g(j 



U-4r-3r0 r 7 r 3 461 r 5 4619r 7 



(A. 2) ib(r) = — I —, = — I 1 1 h 

v 1 vy > r i_ /4-9r2\ 6 40 2016 13440 



log 



From this numerical evidence (see Figure [3] for a plot) it is clear that, for 
r G [0,0.4], 

r 

H r ) ^ and 
D 

V>"(r) 189 r 12862 r 3 44487 r 5 346403009 r 7 189 

+ 77T77^ + . . . < — r. 



(^'(r)) 2 5 175 500 4042500 ~ 5 

Hence £T'(r) = G'(ip(r))ip'(r) < G'(§)^'( r )- As 

= (G o = o ^ + • (V.') 2 • (C o ^ , 

H"(r) ^ follows provided ^ — ^pr. Therefore, assumption (|2. 12 j) 

is fulfilled, if 

For G(x) = j4tanh(^x), in which case G'(x) = B / cosh(-^x) 2 and = 

— 2-j tanh(-^x) this can be checked numerically. (Observe that A ^ 0.4 
and distinguish the cases B = G'(0) ^ 6 and B > 6.) For an illustration see 
the rightmost plot of H(r) in Figure EJ 
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